Search for neutrinoless r decays involving the Kg meson 
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We have searched for lepton flavor violating decays of the r lepton with at least one Kg meson in the final 
state. The data used in the search were collected with the CLEO II and II. V detectors at the Cornell Electron 
Storage Ring (CESR) and correspond to an integrated luminosity of 13.9 fb -1 at the T(4S) resonance. No 
evidence for signals were found, therefore we have set 90% confidence level (C.L.) upper limits on the branching 
fractions B(r~ -> e~ K%) < 9.5 x 1(T 7 , B{r~ -> fj,~ K° 8 ) < 9.5 x 1(T 7 , B(t~ -> e~ K%K° S ) < 2.4 x 1(T 6 , and 
B(t~ — > fi~ K S K S ) < 3.5 x 10 -6 . These represent significantly improved upper limits on the two-body decays 
and first upper limits on the three-body decays. 



1. INTRODUCTION 

In physics, all fundamental conservation laws 
are expected to have associated symmetries. Lep- 
ton flavor conservation, however, is an experi- 
mentally observed phenomena with no associated 
symmetry in the Standard Model. Lepton flavor 
violation (LFV) is expected in many extensions of 
the Standard Model such as lepto-quark, super- 
symmetry, superstring, left-right symmetric mod- 
els, and models that include heavy neutral leptons 
(IJII . Experimentally, both Super Kamiokande || 
and SNO (i) observe neutrino oscillation, which 
may imply LFV in the neutrino sector, therefore 
LFV is expected to occur in charged lepton de- 
cay at some branching fraction, albeit very small. 
The r lepton provides a clean laboratory for such 
searches. Ilakovac || has calculated upper lim- 
its on branching fractions for many neutrinoless 
LFV modes within a model involving heavy Dirac 
neutrinos. The branching fractions depend on the 
heavy neutrino masses and mixings. For the de- 
cays r~ — > t~ Kg ||, the branching fractions are 
of 0(1O -16 ), where I can be e or /i. For the decays 
T~ -► t~K° s K s or t~ -> i~K + K~, the branch- 
ing fractions are of 0(1O~ 7 ). The decays with two 
kaons in the final state are therefore of particular 
experimental interest. Previous published upper 
limits on the branching fractions for the decays 
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t~ -> i~K a s are of 0(l(r 4 ) §. There are no 
previous results for the decays t~ — * £~K S K S . 
In this paper, we present the result of a search 
for the decays into one lepton and one or two Kg 
mesons, with the K s decaying into two charged 
pions. 

The data used in this analysis were collected us- 
ing the CLEO detector [|]|| from e + e _ collisions 
at the Cornell Electron Storage Ring (CESR) at 
a center-of-mass energy y/s ~ 10.6 GeV. The 
total integrated luminosity of the data sample 
is 13.9 fb^ 1 corresponding to the production of 
N TT = 1.27 x 10 7 t+t- events. The CLEO detec- 
tor is a general purpose spectrometer with excel- 
lent charged particle and shower energy detection. 

2. EVENT SELECTION CRITERIA 

The t + t - candidate events must contain four 
or six charged tracks with zero net charge. To 
reject beam-gas events, the distance of closest ap- 
proach of each non-K s track to the IP must be 
within 0.5 cm transverse to the beam and 5 cm 
along the beam direction. The polar angle 9 of 
each track with respect to the beam must sat- 
isfy | cos ^| < 0.90. Photons are defined as en- 
ergy clusters in the calorimeter with at least 60 
MeV in the barrel (|cos0| < 0.80) or 100 MeV 
in the endcap (0.80 < |cos0| < 0.95). We fur- 
ther require every photon to be separated from 
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the projection of any charged track by at least 
30 cm unless its energy is greater than 300 MeV. 
In order to diminish QED background such as 
radiative Bhabha and /z-pair events with photon 
conversion, we require each event to have total en- 
ergy less than 95% of the center-of-mass energy. 
This requirement rejects most of the QED back- 
ground, while incurring a small loss in detection 
efficiency. 

We divide each event into two hemispheres (sig- 
nal and tag), one containing one charged track 
and the other containing three or five charged 
tracks, using the plane perpendicular to the 
thrust axis [|l0). The thrust axis is calculated 
from both charged tracks and photons. The in- 
variant mass of the tag hemisphere must be less 
than the r mass, M T = 1.777 GeV/c 2 [JO). The 
signal hemisphere must contain an electron or a 
muon and one or two K s mesons. The electron 
candidate must have shower energy to momentum 
ratio in the range, 0.85 < E/p < 1.10, and when 
available, the specific ionization lost must be con- 
sistent with an electron. The muon candidate 
must penetrate at least three absorbtion lengths 
of iron. The Kg candidate is reconstructed in 
the 7r + 7r~ final state with a detached vertex, and 
the invariant mass must be within approximately 
three standard deviations of the nominal mass, 
485 < TOtttt < 510 MeV/c 2 . In order to diminish 
radiative Bhabha and /i-pair events further, nei- 
ther pion should be consistent with identification 
as an electron. 

Since there is no neutrino in the signal hemi- 
sphere while there is at least one neutrino unde- 
tected in the tag hemisphere, the missing momen- 
tum of the event must point toward the tag hemi- 
sphere, < cos Qtag-missing < 1-0. In order to 
suppress the background from radiative Bhabha 
and ^-pair events, the direction of the missing 
momentum of the event is required to satisfy 
| cos9 m i SS ing \ < 0.90. For the decay t~ — > e~K s , 
cos Qtag-missing is further restricted to be less 
than 0.99 to reduce the radiative backgrounds. 
This corresponds to the minimum ratio of back- 
ground to detection efficiency. The background 
is estimated from the sidebands in the invariant 
mass vs. total energy distribution of the decay 
candidates. 



To search for decay candidates, we select r can- 
didates with invariant mass and total energy con- 
sistent with the expectations. The following kine- 
matic variables are used to select the candidate 
events: 

AE = E — Ebeam 

AM = M — M T , 

where Eb eam is the beam energy, and E and M 
are the reconstructed r energy and mass. The 
AE vs. AM distributions of the decay candidates 
in the data and signal Monte Carlo samples are 
shown in Figs. [I] and ||. 

In the Monte Carlo simulation, one r lepton 
decays according to two- or three-body phase 
space for the mode of interest, and the other 
r lepton decays generically according to the 
KORALB/TAUOLA r event generator Q. 
The phase space model is appropriate for an un- 
polarized tau. If the Lorentz structure of the 
neutrinoless decay is V-A (Standard Model-like), 
correlations between the spins of the two r's in 
the event will lead to slightly higher detection 
efficiency than phase space, while V+A decays 
will lead to lower detection efficiency. The de- 
tector response is simulated using the GEANT 
program JO . The estimated detection efficien- 
cies (e) |l4| are summarized in Table [|. 



Table 1 

Summary of detection efficiency (with statistical 
error), 90% C.L. upper limits on the branching 
fraction with and without including systematic 



errors. 


Mode 


e(%) 


(stat.) 


B{10-<) 


e-K s 
»-K s 
e-K s K° s 
H-K%Kl 


18.8 ±0.4 
19.0 ±0.4 
12.0±0.1 
8.0 ±0.1 


8.8 
8.7 
21 
31 


9.5 

9.5 

24 
35 
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Figure 1. AE vs. AM distribution of the (a) 
data and (b) signal Monte Carlo sample for the 
decay r~ — > e~K ( g\ (c) and (d) show the cor- 
responding distributions for t~ — ► (i~ Kg. The 
normalization of the signal Monte Carlo sample 
is arbitrary. The ellipses indicate the signal re- 
gion (see text). 



Figure 2. AE vs. AM distribution of the (a) 
data and (b) signal Monte Carlo sample for the 
decay t~ — > e~K s K s ; (c) and (d) show the corre- 
sponding distributions for t~ — ► fi~KgKg. The 
normalization of the signal Monte Carlo sample is 
arbitrary. The ellipses indicate the signal region 
(see text). 
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3. RESULTS 

The upper limit on the branching fraction is 
related to the upper limit A on the number of 
signal events by 



2eN TT B 1 (B K o^ + „-)* ' 

where Bi = (84.71 ± 0.13)% is the inclusive 
1-prong branching fraction Q, B K o^+ n - = 
(68.61 ± 0.28)% is the branching fraction for K s 
to decay to two charged pions, and n is the num- 
ber of Kg mesons in the final state. No candidate 
decays are observed, so we take A as 2.44 events 
in each mode at 90% confidence level, according 
to the frequentist method fig] . The upper lim- 
its on the branching fractions with and without 
systematic errors are shown in Table [l| 

The systematic errors include the uncertainties 
in the r + r~ cross section (1%), luminosity (1%), 
track reconstruction efficiency (1% per charged 
track), Kg detection efficiency (2% per Kg), lep- 
ton identification (1.5% for electron and 4% for 
muon), and the statistical uncertainties in the de- 
tection efficiencies due to limited Monte Carlo 
samples (1-2%). 

Black et al. |ll| have analyzed the con- 
straints on the new physics scale for dimension- 
six fermionic effective operators involving tau-mu 
mixing, motivated by the observed — v T oscil- 
lation. The most stringent lower limits from ex- 
otic heavy quarks and tau decays on the physics 
scale of the operators involving quarks are ~10 
TeV. The new upper limit on B(t~ — > ^ Kg) 
presented in this paper yields a lower limit of 
17.3 and 18.2 TeV for the axial vector and pseu- 
doscalar operators, respectively. 

In conclusion, we have searched for r decays 
involving Kg mesons that violate lepton flavor, 
but find no evidence for a signal. This re- 
sults in improved upper limits for the decays 
t~ — > £~ Kg and first upper- limits for the de- 
cays t~ — > £~KgK s . The upper limits for the 
t~ — > £~KgKg final states are more stringent 
than those found previously for £~K + K~ p7| . 

This work was supported in part by the U.S. 
Department of Energy. 
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